Electron Beam Instability in Left-Handed Media 
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We predict that two electron beams can develop an instability when passing through a slab of left- 
handed media (LHM). This instability, which is inherent only for LHM, originates from the backward 
Cherenkov radiation and results in a self-modulation of the beams and radiation of electromagnetic 
waves. These waves leave the sample via the rear surface of the slab (the beam injection plane) 
and form two shifted bright circles centered at the beams. A simulated spectrum of radiation has 
well-separated lines on top of a broad continuous spectrum, which indicates dynamical chaos in the 
system. The radiation intensity and its spectrum can be controlled either by the beams' current or 
by the distance between the two beams. 

PACS numbers: 41.75.-i, 41.60.Bq 
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Metamaterials (also known as left-handed media 
(LHM) [1]), characterized by having negative both per- 
mittivity e and permeability /i, show numerous re- 
markable counterintuitive features 0, including: nega- 
tive refraction [![, inverted Doppler shift and reversed 
Cherenkov cone (backward Cherenkov radiation) as 
well as surface waves propagating along the interface be- 
tween right- and left-handed media Q ■ These properties 
(promising for applications in a wide frequency range for 
subwavelength imaging and lensing Q , subwavelength Q 
and open Q resonators @], nonradiating configurations 
[§], etc.) have stimulated enormous scientific activity 
during the past decade. However, as far as we know, the 
collective interaction of charged particle beams with elec- 
tromagnetic waves in LHM has not been studied. Nev- 
ertheless, one can expect a very nontrivial interaction 
among particles in the beams since the Cherenkov radi- 
ation emitted by a particle propagates backward 0, [l(| , 
producing strong positive distributed feedback for parti- 
cles moving behind. Such a strong coupling can create 
an instability and chaotic (i.e., very irregular) motion in 
the beams and waves. In this article we predict the insta- 
bility of electron beams in LHM associated with tunable 
self-sustained electromagnetic radiation. 

In order to provide an intuitive picture of the effect dis- 
cussed below, we compare information transport by elec- 
tromagnetic waves (solid line) and particles (dashed line), 
shown in Fig. [1] for slabs of left-handed (a) and right 
handed (b) media. Information about any perturbation 
created, for instance, at the rear (left) surface is trans- 
fered either along the wave characteristics (solid line) 
and/or along the particle trajectories (dashed line) deep 
towards the sample. For right-handed media (RHM), 
both information fluxes are directed forward and per- 
turbations propagate within the shaded region between 
the characteristics (Fig.Q})). Thus, any knowledge about 
the perturbation leaves the sample in the forward direc- 



tion after a finite time. In contrast, for LHM, particles 
and waves transport information in opposite directions. 
Therefore, the information transported from the rear to 
the front surfaces by particles is returned back by the 
emitted waves (Fig. QJi). Simultaneously, these "return- 
ing" waves perturb the particles entering the sample later 
on (in the rear surface), and the process of information 
transport will oscillate back and forth, and never stops. 
Similar processes occur in microwave devices as traveling 
wave tubes (TWT) (similar to a slab of RHM) and back- 
ward wave oscillators (BWO) (similar to a slab of LHM) 



It is well known [ll[ that positive feedback in BWO 
produces an electron beam instability and self-excited mi- 
crowave radiation. Thus, a similar instability can be ex- 
pected in LHM. Consider, for example, a small perturba- 
tion (denoted here as "bunch" ) of the beam density with 
a longitudinal dimension smaller than the wave-length of 
the synchronous wave (i.e., a wave with phase velocity 
equal to the beam particles' velocity). Particles forming 
the bunch radiate coherently and the radiated electro- 
magnetic field is not compensated (due to destructive 
interference) by other beam particles. This wave acts on 
electrons moving behind the bunch and periodically mod- 
ulates their velocities. These velocity variations modu- 
late the beam density resulting in secondary bunches in 
the beam. Due to the periodicity of the emitted waves, 
radiation produced by the bunches sum up coherently 
and the total field increases. A stronger field leads to a 
faster and deeper modulation of the beam density that, 
in turn, strengthens the field. As a result, the amplitude 
of the synchronous wave grows exponentially. This is the 
so-called beam instability, which has been discovered si- 
multaneously in plasma physics [l2l [l3j and microwave 
electronics {14j . Thus, two collective effects (the coherent 
Cherenkov radiation of the particles forming the bunch 
and the coherence between the bunches) constitute the 
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by the standard formula 17 1 



X COS 



— (z — vt) uj du, (2) 
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where v is the particle velocity, (3 = v/c, and n = —yfsji 
is the refraction index. The integration domain in a LHM 
is restricted by the conditions u> > 0, \n\j3 > 1, e < and 
/i < [i.e., < < Up/3/ (1 — /?)]. The other components 
of the radiated fields will not be needed because only the 
i^-component governs the particle motion along the z- 
axis in a strong guiding magnetic field. Below we consider 
the non-relativistic limit, j3 <C 1, and small dissipation, 
v -C f3 1 / 2 oj p . The motion of particles interacting through 
their emitted waves can be described by: 
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FIG. 1: Waves (solid line) and particles (dashed lines) char- 
acteristics in LHM (a) and RHM (b) slabs with thickness L. 
For LHM, particles and waves transport information in the 
opposite direction, while, for RHM, information can only be 
transmitted forward. (c,d) Cherenkov field, emitted by par- 
ticles (dashed lines), located within Cherenkov cones shown 
by solid lines. For case (c) the distance between the beams 
is too large, D > fl max , and the radiation emitted by one 
beam cannot reach the other beam, so there is no interac- 
tion between beams. For case (d), D < Dmax, the radiation 
emitted by one beam reaches the tail of the other beam, cre- 
ating a particle-density modulation producing a strong beam 
coupling. 



physical basis of the beam instability. 

Here we consider two separated beams moving along a 
strong magnetic field and interacting via their Cherenkov 
radiations propagating backward in 3D LHM and predict 
a beam instability. This instability develops when the 
distance between the beams is within a certain interval 
and when the beam currents exceed a certain threshold. 

Model. — We describe LHM by a frequency-dependent 
permittivity and permeability (see, e.g., @, 0, HI, IH): 



e{u>) = fJ,(u>) = 1 



(1) 



where v is the collision frequency and u p is the plasma 
frequency. The medium is "left-handed" when the wave 
frequency lu < u> p . 

We consider two parallel beams, separated by a dis- 
tance D propagating along the z-direction (the limit 
D — > 0, corresponding to one beam, is discussed below). 
Charged particles in these beams emit Cherenkov radia- 
tion with the z-component of the electric field described 
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is the coordinate of the i-th particle in the k- 



k = 1, 2, and E. 
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tj r ), rij] is the radiated field of the j-th particle from the 
fc-th beam acting on the i-th particle from the £-th beam, 
z^} and tj^} are the position and the time when the j-th 
particle radiated the wave that reaches the i-th particle 
at position Zi at time t, ry = when k = £ and = d 
otherwise. E z q(z — f3ct, r) is the "elementary" field with 
space-time structure defined by Eq. (|2|) and depicted in 
Fig. UK. 

The essential property of the Cherenkov radiation 
emitted by a single particle radiation is that the field 
is mostly concentrated along the line p — £ tan a and the 
angle a = arcsin(l/3) is independent of the particle ve- 
locity; in full analogy with the wake that a moving ship 
produces on the surface of the sea 18| . This means that 
we can associate a unique group velocity v g to the largest 
fraction of the radiation. Thus, we can approximate the 

as: 



position Zj r and the time tj r 



tj r — tj - 



d?S 



Vj{ Z 



— , Z jr = Z+\v gz \ (t- tj r ) COS a, (4) 



where tjo is the time when the particle crosses the rear 
surface z — of the LHM slab. 

Hereafter the following dimensionless variables are 
used: £ = k 0z z, r) = k 0z x, ( = k 0z y, and r = k 0z v t 
(vq is the unperturbed beam velocity). 

Instability of two beams. — We simulated the motion 
of electrons in two thin beams, separated by a distance 
D in the transverse direction, using Eqs. ((3]) and {4|; 
each simulation takes time r s 3> £, where £ = k^ z L is 
the dimensionless thickness of the slab. The simulations 
showed that, at certain conditions, a beam instability 
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FIG. 2: (Color online) Spatial distribution of the Cherenkov 
radiation emitted by a single particle, (a) The field distri- 
bution E z in the (p, £')-plane, where p — ko z r and £' = 
koz (z — vt) are the normalized radial and axial distances from 
the particle; ko z is the axial component of the wave vector fco 
corresponding to the maximum of the spectrum of radiation. 
The electric field amplitude is depicted with the factor ^/p, 
which compensates the field decay E ~ p -1 ^ 2 - (b), (c): axial 
distribution of the field E z for p — 12.5 (white line in Fig. 2a) 
and p = accordingly. Note that E z (£) for p = 12.5 has the 
form of a wave packet with a distinguishable periodicity (b), 
whereas E z (£) for p = is fast-decaying and non-periodic (c) . 



develops: any small perturbation develops into a strong 
modulation of the beam density. Inside the slab, the 
beam modulation increases strongly away from the rear 
surface and exhibits a maximum approaching the front 
surface (Fig. [3^)- The beam density modulation is the 
source of radiation which is transported along the re- 
verse Cherenkov cone (see Fig. [5^ and Figs. [TJ:,d) by the 
backward waves to the rear surface. The total field E z is 
calculated as a sum of the "elementary" fields radiated by 
all the particles inside the slab. The total radiation inten- 
sity, hot = J d v d( 1(7], C, £ = 0) = t- 1 J d V d(dT E 2 z (li = 
0, rj, £, r) is not zero when the beam separation A = ko z D 
is within a certain interval A £ (A m ; n , A max ) (Fig. [3Jd). 
The window (A m ; n , A max ) of this instability weakly de- 
pends on the beam current (for simplicity, we consider 
two beams with the same current </{,). The spatial distri- 
bution of the radiation intensity I(rj, £) at the rear surface 
forms two overlapping or intersecting rings, as shown in 
Fig-Et- These bright rings could be seen by an observer 
located outside the slab. 

The physical origin of the suppression of the instabil- 
ity for any large separation distance A > A max is rather 
obvious: indeed, in this case the neighboring beams are 
located outside their corresponding Cherenkov cones and 
cannot interact via their emitted waves (Fig. QJ;). At 
small distances A < A min , the instability is suppressed 
due to three different reasons. First, the electromagnetic 
field emitted by a particle decreases fast (see Fig. [2fc) 
along the beam. Second, the spatial period of the de- 
caying wave varies behind the particle, which destroys 
the spatial periodicity of the beam modulation and, as 



a result, suppresses the coherence of the radiation emit- 
ted by this modulation. These two reasons are common 
for the destruction of the beam instability in both 3D 
LHM and BWO and explain why there is no instability in 
the one-beam system. The other, third reason, is unique 
for 3D LHM: when the distance between the beams de- 
creases, the region of intersection of one beam with the 
Cherenkov cone formed by the second beam is shifted to- 
ward the front surface of the slab. The beam modulation 
is amplified along the distance between this intersection 
and the front surface of the slab. This distance is short- 
ened when the beam separation decreases and becomes 
insufficiently long for an instability to develop. Particles 
in the two different beams interact effectively and the ra- 
diation intensity reaches its maximum when the intersec- 
tion between the beam and the Cherenkov cone occurs 
near the front surface (Fig. QJ1). This is in agreement 
with our simulations (Fig. [3^,). Therefore the distance 
between the beams can effectively control the intensity 
of the radiation. 
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FIG. 3: (Color online) (a) Amplitude of the beam density 
modulation (due to bunching of particles) along the beam, 
divided by the average density. The maximum of the density 
modulation occurs near the front surface; thus, beams radiate 
with higher intensity near this surface, (b) Total radiation 
intensity 7 to t versus the separation distance A between two 
beams. This intensity is nonzero for A m i n < A < A max . (c) 
Spatial distribution of the radiation intensity I(rj, f ) at the 
rear surface of the LHM slab. This distribution has the form 
of two shifted circles centered at the beams. 



Properties of radiation. - Let us consider the ra- 
diation for the optimal separation between the beams. 
An analogy with BWO [ll| points out that we can ex- 
pect three different radiation regimes: (i) no radiation 
for either a weak beam current or small slab thick- 
ness, (ii) monochromatic radiation above a certain cur- 
rent/thickness thresholds, and (iii) periodic or chaotic 
self- modulation at larger currents/thickness. Indeed, we 
found no instability at either small current or thin sam- 
ples. Above a threshold we always observed a chaotic 
self-modulation. An example of the time-dependence of 
the electric field of the emitted wave at the rear surface is 
presented in Fig. [JJi. This dependence clearly shows the 
irregular modulation of the wave amplitude. The spectra 
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of the wave is shown in Fig. [4] (b-d) , for increasing beam 
current. All of these three presented spectra have a simi- 
lar structure: narrow spectral lines whose overlap forms a 
continuous spectrum. With increasing current, the lines 
become broader and the strength of the continuous spec- 
trum grows, which is usual for microwave devices with 
a so-called "frequency" scenario of transition to chaos 
[Tit . The frequency gap between the spectral lines is de- 
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FIG. 4: (a) Time-dependence of the electric field E(t) of the 
emitted wave at the rear surface of the LHM slab. The irreg- 
ular modulation of the wave amplitude indicates chaotic dy- 
namics of the beams, (b-d) Radiation spectra of a two-beam 
system when the beam current grows. Spectra show well sep- 
arated narrow lines on top of a continuous background. The 
current increases by a factor 1.5 from (b) to (d). 

termined by the time needed for a signal to close the 
feedback loop (signal is transported first by one beam, 
then the signal is carried back by the electromagnetic 
wave from the beam head to the tail of the second beam, 
then forward along the second beam, and afterward by 
the wave from the head of the second beam to the ini- 
tial point; see inset in Fig. UJd). Therefore, the distance 
between the spectral lines can be tuned by the distance 



between the beams or can be designed by a proper choice 
of the slab thickness. 

Conclusions. — We predict a beam instability in 3D 
left-handed media. This instability produces strong ra- 
diation which can be tuned either by the beam current 
or by geometrical parameters. The spectrum of the ra- 
diation has well separated lines on top of a broad back- 
ground. Since now LHM can be fabricated [2(| in the 
optical frequency range, the predicted effect can be used 
to generate and amplify optical irregular signals. Further 
extension of the analogy between traditional microwave 
devices and left-handed media will allow new insights for 
new potential applications of LHM. 
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